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&-*'■ Abstract 

We present candidates for the global minimum energy solitons of charge one to nine 
in the Skyrme model, generated using sophisticated numerical algorithms. Assuming 
the Skyrme model accurately represents the low energy limit of QCD, these config- 
urations correspond to the classical nuclear ground states of the light elements. The 
solitons found are particularly symmetric, for example, the charge seven skyrmion 
has icosahedral symmetry, and the shapes are shown to fit a remarkable sequence 
defined by a geometric energy minimization (GEM) rule. We also calculate the en- 
ergies and sizes to within at least a few percent accuracy. These calculations provide 
the basis for a future investigation of the low energy vibrational modes of skyrmions 
and hence the possibility of testing the Skyrme model against experiment. 



The Skyrme model presents an opportunity to understand nuclear physics as a low 
energy limit of quantum chromodynamics (QCD). It was first proposed as a theory of 



strong interactions of hadrons[16], but after a hiatus which lasted nearly twenty years it 



was shown to be in fact the low energy limit of QCD in the large N c limit |17|. Since then 
much work has suggested that the topological solitons in the model, known as skyrmions, 
can indeed be identified with classical ground states of light nuclei. However, before a more 
quantitative assessment of the validity of the model can be made, a thorough understanding 
of the structure and dynamics of multi-soliton configurations is required. 

This is a highly non-trivial task since the Skyrme model has many features which ob- 
struct both analytic and numerical progress. There are two promising lines of attack. 
The first is to use the approximate correspondence between Skyrme fields and SU(2) 
instantons[|l|]. Using this semi-analytic approach, one can understand the right-angle scat- 
tering of two skyrmions0 and also calculate symmetric configurations of, for example, 
charge three and fourjnj. The second is to use state of the art computer technology and 



a sophisticated numerical algorithm to evolve the dynamics of discretized configurations. 
This approach has the advantage of being quantitatively accurate, but has many other 
more practical difficulties associated with it, in particular the substantial computer re- 
sources required. Nonetheless, it has been possible to design a working code which can 
evolve the dynamics of low energy configurations for many timesteps and also be modified 
to relax dynamically towards what are at least local minimum energy configurations. 

In a recent letter |J, we presented the first results from our code, exhibiting the dynamics 
of symmetric configurations in the attractive channel and also the minima for charges one to 
four. These minima were already known from numerical relaxation calculations ]7| and also 
from instanton calculations [[TJ]]. In this letter, we present what we describe as candidate 
minima for charges one to nine, which we firmly believe are in fact the global minima. These 
minima can be classified in terms of the isosurfaces of baryon density (or energy density), 
which are seen to fit a remarkable pattern, obeying what we shall call a geometric energy 
minimization (GEM) rule. We also calculate accurately the energy and average size of 
each soliton. As in ref. ||, we shall assume the pion mass is zero -- results for a physical 
value of the pion mass will be presented in future work. We do not expect that inclusion 
of a finite pion mass will effect the overall shape and symmetry of the soliton, rather it will 
just modify its energy and size. 

The Lagrangian of the massless Skyrme model may be written in terms of the SU(2) 
valued right currents R^ = (d^U)W as 

12tc 2 £ = -±Tr(iV^) " ±Tr{[R ll ,R l/ ][R»,K']) 

where U = a + it ■ tt and we have used scaled units of energy, E = 661MeV, and length, 
£ = 0.755fm. One can convert to an 0(4) sigma model representation <p = (a, 7i"i, 1T2, K3) as 
described in refs. ||, |6| and deduce the dynamical equations of motion for <fi. The baryon 
density B, whose spatial integral gives the integer- valued baryon number B, is given by 
2Att 2 B = —eijkTi{RiRjRk) where Latin indices run over the spatial values 1, 2, 3. The 
above units are chosen so that the Fadeev-Bogomolny bound on the energy E is simply 
E> \B\. 



The full details of the numerical methods are presented in reffl. This includes an 
account of how the numerical algorithm works, how to construct the initial conditions and 
why we believe that the configurations represent the global minima. For the purposes 
of this letter, we include a brief outline for the layman. The basic procedure is to set 
up discretized initial conditions for static low energy configurations using either a simple 
ansatz for a single skyrmion and the product ansatz, or by calculating the holonomy of 
instantons, or a combination of both, as is the case for the higher charge configurations. 
The critical feature of the specific configurations chosen is that the skyrmions are in an 
attractive or nearly attractive channel. These configurations are then evolved using the 
discretized equations of motion. As the dynamics proceeds the system oscillates between 
maxima and minima of the potential energy. If one stops the dynamics at a minimum 
of the potential energy and then removes all the kinetic energy from the system, that is, 
set the time derivatives of the field to zero, the solution will gradually move towards a 
local minimum of the system. Once the system is close to the minimum one can also 
incorporate dissipation which will speed up the process of relaxation. Of course, as with 
any numerical minimization procedure, one cannot be sure that one has found the global 
minimum. However, given sufficiently asymmetric, but attractive, initial conditions, as 
created using the product ansatz, it is likely that one will locate the global minimum. 

For charges one to four it was possible to construct configurations using just the product 
ansatz, in which all the skyrmions are in a mutually attractive channel^, |6|. It is not 
possible to use this naive approach to construct a maximally attractive channel of charge 
five or higher and therefore one must consider other approaches. In order for our algorithm 
to relax quickly to the minimum, we must construct a configuration with low energy, in 
which most of the skyrmions are attractive. It need not be the maximally attractive 
channel, but the algorithm will relax to the minimum quicker if it is close to the most 
attractive. Some of the configurations chosen relax quicker than others, suggesting that 
in some cases we have selected the correct configuration and in others we have found one 
which works, but is perhaps not the best. 

We find that it is best to mix the instanton approach for a known symmetric configura- 
tion, with a small number of single skyrmions added using the product ansatz to break the 
exact symmetries. For example, in the case of B=5, it possible to construct a highly at- 
tractive configuration by adding in two single skyrmions either side of a B=3 tetrahedron. 
Since the pion fields of the tetrahedron are similar to those of an anti-skyrmion at large dis- 
tances, the configuration is of low energy and relaxes extremely quickly (~1000 timesteps) 
to the minimum. Using the analogy, to B=2 to B=4 scattering, we tried constructing 
B=6 and B=7 configurations by surrounding a B=3 tetrahedron with skyrmions in cyclic 
and tetrahedral configurations respectively using the product ansatz. These configurations 
relaxed toward the minimum in over 5000 timesteps, suggesting that in fact these are not 
the maximally attractive channel. Nonetheless, they are of reasonably low energy and work 
eventually. For B=Q, one can relax to the minimum much more quickly by colliding two 
tetrahedra. Reassuringly, it is the same minima as calculated before, providing a useful 
consistency check. In order to calculate configurations of higher charge (B > 7), one can 
just add in a single skyrmion to the known minima of charge one less. 



A useful way to represent a skyrmion is by displaying a surface of constant baryon 
density. In fig. 1 we display isosurfaces for the skyrmions of charge five to nine (one to four 
are presented in ref. 0), using the same constant value for the baryon density in each case, 
to ensure that the relative sizes of the skyrmions are accurately represented. The baryon 
density has maxima at several points in space, which we can think of as vertices of a solid, 
and these are connected by links of slightly lower baryon density, which can be regarded 
as the edges of the solid. Clearly in this way we can assign a solid to each skyrmion which 
accurately reflects the shape and symmetry of the skyrmion. In fact, for all the skyrmions 
we consider, the associated solid is remarkably close to being composed of regular or nearly 
regular polygons with a fixed edge length. Given this fascinating result, we shall describe 
each solid in terms of its construction from polygons. Included in fig. 1 alongside each 
baryon density plot is a photograph of the associated solid, constructed using a molecular 
model building kit. In table 1 we list the number of faces of each associated solid, together 
with the number of each type of polygon from which it is constructed. 
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Table 1 : The constituent polygons and symmetry structure for the candidate minima 
of charge three and above. 

It is perhaps useful to give a brief description of each of the solids, which in conjunc- 
tion with fig. 1, should make the structure of each configuration clear. The B=5 skyrmion 
consists of two parallel down-pointing pentagons attached to two more parallel up-pointing 
pentagons, so that four sides of a box are formed. The top of the box is formed by adding 
two squares, and similarly for the bottom of the box, though of course the arrangement of 
joined squares on the top and bottom have a relative rotation of 90°. The B=6 configura- 
tion consists of two halves, each of which is formed from a square with a pentagon hanging 
down from each of its four sides. Note that to join these two halves implies that the two 
squares are parallel, but one is rotated by 45° relative to the other. The B=7 solid is a 
regular dodecahedron. The B=8 skyrmion has a similar structure to its B=6 counterpart, 
except that the squares are replaced by hexagons, so that each half has six pentagons 
hanging down. This requires the top hexagon to be parallel to the bottom hexagon but 
rotated by 30°. Finally, the B=9 structure has four hexagons located at the vertices of 
a regular tetrahedron which are joined by four sets of three connected pentagons, whose 
single common vertex lies at the vertices of the dual tetrahedron. 

The symmetries of each skyrmion are of interest and here we shall discuss the symmetry 
group of the baryon density (or equivalently energy density) . To discuss the symmetries of 



the field itself is a more complicated task, since the three-dimensional representation of the 
symmetry group which acts on the pion fields as an isospin rotation must also be identified. 
In the final column of table 1 we give the symmetry group, where we use the Schonflies 
notation (see for example, ref.[15|) popular in chemistry. All the configurations contain 



at least the symmetry group D n d : for some n. This symmetry group is obtained from the 
cyclic group of order n, C n , by the addition of a Ci symmetry with axis perpendicular to 
the main symmetry axis and a reflection symmetry in a vertical plane containing the main 
axis and which bisects the angle between pairs of C% axes. Thus these twisted dihedral 
symmetries (which are enhanced to platonic symmetry in some cases) appear to be of 
importance to skyrmions, as they are for BPS monopoles |[T0|, |TT|| . 

Note that the skyrmions of charge 3,4,7 have the same platonic forms as the corre- 
sponding BPS monopoles of the same charge[RL pi fLO], but that the charge 5 configuration 



is not an octahedron, even though an octahedral 5-monopole exists[10|. An octahedral 
charge 5 skyrmion does exists but it is not the minimal energy configuration|B[]. This illus- 
trates an important difference between BPS monopoles and skyrmions. All the monopole 
configurations of a particular charge have the same energy and it is due to a mathematical 
simplification that only the very symmetric ones have been found, but the skyrmions evolve 
under the influence of a potential allowing for the possibility of minima. Nonetheless, it is 
fascinating that in reality the highly symmetric low energy configurations of B=3, 5=4, 
B=7 and B=9 form in the early universe during Big Bang Nucleosynthesis, while the less 
symmetric configurations of B=5, B=6 and B=8 are either unstable or barely stable. 

Given the complex structure of these skyrmions the question now arises as whether a 
rule exists which fits the remarkable sequence of shapes found. We propose the following 
phenomenological rule for the structure of the minimum energy charge B > 2 skyrmion, 
which we refer to as the Geometric Energy Minimization (GEM) rule; 

GEM RULE : The charge B baryon density surface is composed of almost regular polygons 
and consists of 4(B — 2) trivalent vertices. If more than one such solid exists, then select 
the most spherical. 

It should be noted that there are several equivalent ways in which the GEM rule could 
be stated. For example, from the trivalent property together with Eulers formula, fixing 
one of the three parameters of the solid, that is, the number of vertices, faces and edges, 
determines the other two. Explicitly we have V = A(B - 2), F = 2(B -1),E = 6(B - 2). 
Thus since the baryon density isosurface has a hole in the centre of each face then the 
GEM rule implies the observation of ref.[[jj that the isosurface contains 1{B — 1) holes. As 
the value of B increases the number of possible configurations satisfying the first part of 
the GEM rule grows, and hence the need for the second part. For example, at B = 6 in 
addition to the configuration found there is a second possibility that consists of a hexagon 
with alternating squares and pentagons rising from its edges, and topped by three joined 
pentagons. This crown-like configuration, which has cyclic C3 symmetry, satisfies the first 
part of the GEM rule, but is not very spherical, having a flat bottom and a pointed top. 
It may be that some other statement, such as a minimization of the standard deviation 



of edge lengths, or a similar property for edge angles, is an improved statement of the 
spherical property. More examples of higher charge skyrmions are required to resolve this 
issue. 

It follows from the GEM rule that for B > 7 the solid consists of 12 pentagons and 
2(5 — 7) hexagons. Such configurations occur in fullerene chemistry P, [I2|| , where we should 
compare with the carbon structure C^b-2) ■ in fullerine chemistry avoiding large curvature 
is important, so the first fullerene is C 2 s, avoiding a fused quartet of pentagons, which 
has precisely the form of the B=9 skyrmion. It would be interesting to see if this corre- 
spondence continues, since the 5=17 skyrmion should have the Ceo Buckminsterfullerene 
structure. Extension of our current results will require considerable amounts of computer 
resources, but this is now in progress. 
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Table 2 : Calculated values of the charge (5 d i S ), energy (-Edis) and soliton size (Ar) in 
natural and physical units for B = 1 to B = 9. Also presented is the ionization energy /, 
which is the energy required to remove one skyrmion. For a discussion of the accuracy of 
the results and comparison to previous calculations see the text. 

We have calculated the discrete charge and energy for the minima on 100 3 grids, which 
are displayed in table 2. These energy values are less than the true values since the grid 
size is finite, but they are nonetheless within 2% accurate. One can make a better estimate 
of the overall energy by using the ratio Edi S /Bdi S as discussed in ref.JBj. This can be seen 
to be exact to 3 decimal places for the B=\ skyrmion and we suggest that it will be so for 
all the othersQ since the third decimal place has not changed for many (>1000) timesteps. 
These values have been used to calculate the ionization energy (7g = Ei+Eb-i — Eb), that 
is, the energy required to remove one skyrmion, which in general gives an indication of the 
classical binding energy. Taking into account the fact that these energies will be modified 
by the quantisation of zero modes and also some low energy vibration modes, there is a 
remarkable trend which follows the pattern of light element binding energies. We have 
also calculated the average size of the soliton Ar from the second moment of the baryon 



1 It should be noted that the quoted values for B=l to B=4 differ very slightly from those presented in 
ref. B. The current values are the result of further relaxation of the same configurations. 



distribution as in ref.0. This value is extremely rough since it ignores the symmetries of 
the object, but nonetheless it gives an indication of the overall trend. 

We should comment on the relevance of our work to previous calculations. In numer- 
ical work0, similar to ours, configurations up to charge six were studied using a global 
minimization algorithm. Our results for charges five and six differ from these earlier compu- 
tations, which we attribute to numerical effects in the earlier work due to lack of resources; 
for example, the grid we use contains almost five-times the number of points. For charge 
five the difference is small, we obtain the same symmetry group but identify different poly- 
gons forming the solid, which is possible thanks to our improved grid resolution. For charge 
six our results are very different, as we find a very symmetric configuration, whereas the 
earlier computation gave a structure with very little symmetry. 

A different approach to constructing high charge configurations has used the Skyrme 
crystal 1 1^1 and involves cutting out sections @]. Although these configurations have low 



energy, they are not as low as those presented here, and they are fundamentally different 
in nature. The configurations presented here are shells with less baryon density in the 
centre, whereas those created from the Skyrme crystal have internal structure since they 
are created from cubic configurations. Such structures do not fit the GEM rules since they 
are not trivalent. It is an open question as to whether the shell structure persists for higher 
charge. 

We believe that the candidate minima which we have presented here are in fact the 
global minima since the initial conditions have natural asymmetry, we have in some cases 
got the same configuration from two different initial conditions and most of all the isosur- 
faces of the baryon density fit a remarkable sequence of polygons. The symmetry proper- 
ties of these polygons could be the starting point for an understanding of the moduli space 
structure of the Skyrme model, leading to a study of low energy dynamics. The physical 
properties of these classical nuclear ground states follow the trend of light elements, at 
least qualitatively, providing the impetus for a study of the low energy vibrational modes 
and the potential for a confrontation between the Skyrme model and nuclear experiments. 
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Baxter, Paul Shellard, Jonathan Moore, Dick Hughes- Jones, Neil Turok and Kim Baskerville. 
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Figure Captions 

Fig. 1. skyrmions of charge five to nine; on the left baryon density isosurfaces (to scale) 
with five at the top and nine at the bottom and on the right wireframe models of the 
corresponding solids (not to scale). 
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